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ABSTRACT 

Context. Identification of y-ray-emitting Galactic sources is a long-standing problem in astrophysics. One such source, 1AGL J2022+4032, 
coincident with the interior of the radio shell of the supernova remnant Gamma Cygni (SNR G78.2+2.1) in the Cygnus Region, has recently been 
identified by Fermi as a y-ray pulsar, LAT PSR J2021 +4026. 

Aims. We present long-term observations of 1AGL J2022+4032 with the AGILE y-ray telescope, measuring its flux and light curve. 
Methods. We compare the light curve of 1AGL J2022+4032 with that of 1AGL J2021+3652 (PSR J2021+3651), showing that the flux variability 
of 1AGL J2022+4032 appears to be greater than the level predicted from statistical and systematic effects and producing detailed simulations to 
estimate the probability of the apparent observed variability. 

Results. We evaluate the possibility that the y-ray emission may be due to the superposition of two or more point sources, some of which may be 
variable, considering a number of possible counterparts. 

Conclusions. We consider the possibility of a nearby X-ray quiet microquasar contributing to the flux of 1AGL J2022+4032 to be more likely 
than the hypotheses of a background blazar or intrinsic y-ray variabilty of LAT PSR J202 1+4026. 



O Key words. Gamma rays: stars: 1AGL J2022+4032, 1FGL J2021.5+4026, 3EG J2020+4017 — pulsars: LAT PSR J2021+4026 
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1. Introduction 

Identification of the Galactic sources emitting y-rays with ener- 
gies above 100 MeV is a long-standing problem in astroph ysics. 
With the laun ch of the AGILE y- ray telescope in 2007 ([Tkvani 



et al. 



2009a) and Fermi in 2008 ( |Atwood et aLj [2009] ), a great 
deal of progress has been made. In particular, Fermi and AGILE 
have shown that the vast majority of Galac tic y-ray sources ar e 
probably pulsars and pulsar wind nebulae (Ab do et aL| [2009a). 
However, new source classes have also been identified, including 
micr oquasars such as Cygnus X-l ( Sabatini et al. , 2010 ) , Cygnus 
X- 3 ([Tavani et al.| |2009ct |Abdo et al.||2009c), LS5039 (|Abdo et| 
|aL| |2009d| ) andLSI+61°303 (LSI+61 303 ) (|Abdo et al.| |2009e), 
and Wolf-Rayet stars such as Eta Carinae (Tavani et al. , 2009b). 

As the identification of the remaining unidentified EGRET 
sources and the newly found AGILE and Fermi-LAT sources 
continues, the problem of source confusion will become increas- 
ingly serious. Although the angular resolution of both AGILE 
and Fermi are higher than that of EGRET, they also are sensitive 
to much lower fluxes, Fermi-LAT in particular. As a result, y-ray 
error contours will not only continue to contain many plausible 
source counterparts, but will be increasingly likely to contain 
multiple real y-ray sources above the nominal flux threshold. 
In this paper we provide evidence that 1AGL J2022+4032 



i i J2. rc 

in the Cygnus region y\{ Pittori et al. 2009| may be one of these 




1 Throughout this paper we refer to the Cygnus region, although 
physically this direction lies along the tangents of at least three different 



Fig. 1. SNR Gamma Cygni (G78.2+2.1) in galactic coordinates. 
DRAO Radio telescope, wavelength=21.1 cm. White contour 
levels: AGILE-GRID intensity contour levels - related to Figure|2] 
- (pixel size 0.1°), starting from 0.00085 in steps of 0.00002 (in- 
tensity per pixel); green contour: AGILE-GRID 95% confidence 
level for E > 100 MeV; black circle: LAT PSR J2021+4026. 
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sources. We present a study of the time variability of the y-ray 
flux of 1AGL J2022+4032 based on long term observations of 
the source by AGILE. This source was first discovered by COS- 
ES as 2CG078+2, and was listed as the unidenti fied source 3EG 
J2020+4017 in the 3rd EGRET catalog ( [Hartman et aT][T999] ). 
It is located (Figure [TJ within the wide radio shell of the promi- 
nent SNR Gamma Cygni (SNR G78.2+2.1). Searches for X-ray 
counterparts of the source were performed with a number of X- 



3. Discussion 

3.1. Variability 



ray telescopes, including ROSAT ([Brazier et al. 1996), ASCA 



( [Uchiyama etaL] [20 02 b INT EGRAL ([Bykov et al.| |2004[ ), and 



Chandra ( [Becker et al.[|2004^|Weisskopf et al.[|2006| ). Using a 
blind search technique, Fermi-LAT was finally able to identify 
this source, 1FGL J2021.5+4026 ( |Abdo et al.[ [2009a] ~ 



2010a) 



2009b). 



as y-ray pulsar LAT PSR J202 1+4026 ( [Abdo et al.f 
The y-ray flux and spectra are consistent with a pulsar origin. 
However, y-ray pulsars show no variability over time scales of 
weeks to months; indeed, significant variability on these time 
scales would be difficult to reconcile with theoretical models. 

Beginning in November 2007, AGILE data from this source 
have reveale d signs of variable y-ray emi ssion, as reported in 
ATel #149 2 ([Longo et al.| |2008] ), #1547 ( [Giuliani et alT|[2QQ8] ) 



and #1585 ( [Chen et al.[|20()8| ). The evidence for y-ray variability 
combined with its relatively high unpulsed y-ray fraction could 
indicate either the presence of an additional y-ray source coinci- 
dent with 1FGL J2021. 5+4026 or variability of the y-ray pulsar 
LAT PSR J202 1+4026. 



2. Observations and Data Analysis 

AGILE has been in orbit since April 2007 and has observed the 
Cygnus region numerous times, beginning in November 2007. 
We used these observations to characterize the y-ray variabil- 
ity of the sources in the Cygnus region using the following 
procedure. First, we used the software tool AG_multi2, one of 
the AGILE Scientific Tools, to perform multi- source likelihood 
analysis on the deep-integration AGILE-GRID data using the 
FM3. 119_2 filter. This analysis revealed y-ray emission from 
four point sources: 1AGL J2022+4032, 1AGL J2021+3652, 
1AGL J2032+4102 and a persistent faint source consistent 
with the position of Cygnus X-3 (Figure [2]). The average 
fluxes and positions are shown in Table [lj Fermi source 
1FGL J2020.0+4049 was not detected in either energy range due 
to the combination of its proximity to 1 AGL J2022+4032 and its 
low flux. 

Next, we divided the observations from November 2007 
to August 2009 into 42 discrete fixed-length time intervals of 
« 6 days (« 90 orbits) each, analyzing the AGILE y-ray flux 
from the position of 1AGL J2022+4032, (l,b) = (78.23, 2.12), 
for E > 100 MeV and E > 400 MeV, while keeping its po- 
sition fixed and all nearby sources fixed in flux and position. 
We performed the same analysis on the nearby y-ray source 
1AGL J2021+3652 3.5° away from 1AGL J2022+4032, which 
|Halpern et alJpOQg) identified as PSR J2021+365, in order to 
account for the effects of systematic errors. Figure [3] shows the 
light curves of the two sources. 

For E > 400 MeV, only the emission from the three steady 
sources shown in the lower half of Table [T] were considered; 
J2033+4050 is not significantly detected for E > 400 MeV. In 
Figure[4]we compare the light curves of 1 AGL J2022+4032 and 
1AGL J2021+3652. 



We used the method developed by McLaughl in et al.| ( |1996| ) 
to test the y-ray flux variability of 1AGL J2022+4032 with re- 
spect to 1AGL J2021+36 52. A similar analysis fo r all the 1AGL 
sources is in preparation ( Ver recchia et al.|[2010| ). The weighted 
mean flux is calculated from the fluxes in each 6-day time inter- 
val and their corresponding errors, from which the^ 2 is derived. 
Q is the probability that an intrinsically non- variable source (i.e. 
with constant flux) would produce by random chance a mea- 
sured value of x 2 greater than or equal to the x 2 observed, and 
the variability index V is defined as V = - log Q. A source can 
be classified as nonvariable if V < 0.5, uncertain if 0.5 < V < 1, 
or variable if V > 1 . The value V = 1 corresponds to a prob- 
ability of variability P var = 1 - Q of 90%. Table [2] shows the 
value of V for 1AGL J2022+4032 and 1AGL J2021+3652 for 
E > 100 MeV and E > 400 MeV. For 1AGL J2022+4032 we 
find V = 2.18 when systematic effects are included (V = 3.88 
for statistical only). For 1AGL J202 1+3652 the corresponding 
values are V = 0.30 (V = 0.48). 

As a cross-check, we also calculated a complementary vari- 
ability ind ex, W, according to the formula used in the Fermi 



catalogs ([Abdo et aL 



2009a 



Galactic spiral arms at different distances, from the local « 1 kpc Orion 
spur, to the Perseus arm, to the Outer arm at > 8 kpc. 



2010a[ ). This index is a simple x 
where the weights include the systematic error, f re i which in our 
case is 10%, and the number of degrees of freedom is 41. The 
values, associated probabilities, and fractional excess variabil- 
ities above systematic and statistical fluctuations are found in 
Tabled 

We find evidence for variability for E > 100 MeV in the 
emission from 1AGL J2022+4032 even allowing for systematic 
errors on the level of 10%. Any systematic effects that would in- 
fluence the measurement of the flux should also have affected the 
nearby source 1AGL J2021+3652, for which no corresponding 
variability is found. However, 1AGL J2021+3652 is only half 
as bright as 1AGL J2022+4032. Similarly, although we found 
no evidence for variability in the flux of 1AGL J2022+4032 for 
E > 400 MeV when systematic errors are taken into account,the 
average flux is only a quarter that of E > 100 MeV. In both cases, 
the same intrinsic variability might be rendered undetectable be- 
cause of reduced photon statistics. 

To determine whether the variability of 1AGL J2022+4032 
for E > 100 MeV would have been detectable if the source 
were half as bright (Case 1) or at higher energies (Case 2), we 
produced simulated observations of hypothetical sources with 
the same intrinsic variability, where "intrinsic" includes the ef- 
fects of systematic errors. The meaning of "same" is some- 
what ill-defined, since many qualitatively and quantitatively dif- 
ferent parent distributions can produce the same variability in- 
dex, or even the same observed fluxes within errors. In order 
to simplify the relationship with the calculation of x 2 \ we as- 
sume that the intrinsic flux in each observation is drawn from a 
Gaussian distribution whose mean is equal to the mean flux of 
the source in order to find the variance which reproduces the 
observed value of V = 3.88 for 1AGL J2022+4032. We de- 
termined that, given the exposures of each of the 42 time in- 
tervals, an intrinsic variability of 26% (square root of variance 
33.8 x 10" 8 photons cm" 2 s" 1 for 131 x 10" 8 photons cm" 2 s"\ 
the mean flux of 1AGL J2022+4032 for E > 100 MeV) pro- 
duces a distribution in V with a median equal to the observed 
value, 3.88 (Figure [5] solid black). 

For Case 1, we simulated 10000 series of 42 observations 
for a source with intrinsic fluxes taken from a Gaussian distribu- 
tion with mean equal to 60 x 10" 8 photons cm" 2 s" 1 and square 
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Fig. 2. Cygnus Region in galactic coordinates, y-ray intensity map for E > 100 MeV. Deep integration AGILE-GRID data 
(November 2007 - August 2009). Pixel size = 0.1° with 3-pixel Gaussian smoothing. Green contours: AGILE-GRID 95% confidence 
level; Cyan contours: AGILE-GRID statistical + systematic error (0.1°); Black contours: Eermi-LAT (1-year catalog), statistical error 
only, crosses shown for contours too small to be visible. The Green contours have been calculated with a multi-source likelihood 
analysis, using four persistent sources (Table [I]). 

Table 1. y-ray sources in the Cygnus Region. 



Name 




Position 






Vts 


Flux a 


1AGL J2021+3652, E > 100 MeV 


0,b) = 


= (75.22, 0.24) 


+ 0.08° 


(stat) 


+ 0.10° (syst) 


25.02 


60 + 3 (stat) + 10% (syst) 


1AGL J2022+4032, E > 100 MeV 


a,b> = 


= (78.23,2.12) 


+ 0.06° 


(stat) 


+ 0.10° (syst) 


39.64 


131 + 4 (stat) + 10% (syst) 


1AGL J2032+4102, E > 100 MeV 


a,b> = 


= (80.08, 1.18) 


±0.18° 


(stat) 


+ 0.10° (syst) 


10.82 


37 + 4 (stat) + 10% (syst) 


J2033+4050^ , E > 100 MeV 


a, b) = 


= (79.84, 0.50) 


+ 0.29° 


(stat) 


+ 0.10° (syst) 


5.17 


14 + 3 (stat) + 10% (syst) 


1AGL J2021+3652, E > 400 MeV 


a,b> = 


= (75.16, 0.23) 


+ 0.07° 


(stat) 


+ 0.10° (syst) 


23.07 


17 + 1 (stat) + 10% (syst) 


1 AGL J2022+4032, E > 400 MeV 


a,b> = 


= (78.21,2.12) 


+ 0.05° 


(stat) 


+ 0.10° (syst) 


33.80 


33 + 1 (stat) + 10% (syst) 


1AGL J2032+4102, E > 400 MeV 


a,b) = 


= (80.05, 0.98) 


+ 0.14° 


(stat) 


+ 0.10° (syst) 


10.59 


9 + 1 (stat) + 10% (syst) 



Notes. y-ray fluxes in units of 10 8 photons cm 



(b) 



Positionally consistent with Cygnus X-3. 



root of variance ofl5.5xl0 8 photons cm 2 s 1 using the same 
exposures as 1AGL J2022+4032 for E > 100 MeV (Figure [5] 
dotted red). The median value of V is 0.50. V < 0.48 is pro- 
duced in 49% of trials, while V > 1.0 in 28% of trials. The same 
level of intrinsic variability would be likely to produce a vari- 
ability index similar to that of 1AGL J2021+3652, but would 
not be likely to be classified as variable, if the source were half 
as bright. For Case 2, we simulated 10000 series of 42 obser- 
vations for a source with intrinsic fluxes taken from a Gaussian 
distribution with mean equal to 33 x 10" 8 photons cm -2 s" 1 and 
square root of variance 8.4 x 10" 8 photons cm -2 s" 1 , using the 
exposures of 1AGL J2022+4032 for E > 400 MeV (Figure [5] 
dashed green). The median value of V is 1.75. V < 1.11 is pro- 
duced in 31% of trials, while V > 1.0 in 74% of trials. The same 
level of flux variability at high energies would be detected more 
often than not. Equivalent intrinsic variability at the flux level 



of 1AGL J2021+3652 for E > 400 MeV would be completely 
undetectable (Figure [5| dot-dashed blue). 



3.2. Analysis of pulsed and unpulsed components 

We attempted to perform an analysis of the on- and off-peak 
components of the 1 AGL J2022+4032 AGILE data with respect 
to the Fermi y-ray ephemeris of LAT PSR J202 1+4026 in order 
to test whether the apparent flux variability is due to a source 
other than the pulsar. However, as shown in Table 3 of the Fermi 
Pulsar Catalog ( |Abdo et al.| |20T0c1 ), LAT PSR J2021+4026 is 
one of only two Fermi pulsars for which the y-ray phase pro- 
file is difficult to separate into on-peak and off-peak phases, both 
because the y-ray peaks of the pulsar are very broad and be- 
cause the y-ray emission has a high unpulsed fraction (Figure[6]). 
We attempted a number of cuts on the AGILE data. Varying the 
off-peak exposure fraction from 10% to 45% yielded a mono- 
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1AGL J2022+4032 (90 orbits data) 
1AGL J2022+4032 (weighted mean flux) 
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MJD 



1AGL J2021+3652 (90 orbits data) 
1AGL J2021+3652 (weighted mean flux) 




54400 54500 54600 54700 54800 54900 55000 55100 
MJD 



Fig. 3. Flux for E > 100 MeV on six-day time intervals from 1 AGL J2022+4032 (left) and 1 AGL J2021+3652 (right). The Red line 
indicates the weighted mean of the 42 six-day fluxes, from which the^ 2 calculated was calculated: (1 18±4)x 10" 8 photons cm -2 s" 1 
for 1AGL J2022+4032, and (54 ± 3) X 10" 8 photons cm -2 s" 1 for 1AGL J2021+3652. 



1AGL J2022+4032 (90 orbits data) 
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1AGL J2021+3652 (90 orbits data) 
1AGL J2021+3652 (weighted mean flux) 
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Fig. 4. Flux for E > 400 MeV on six-day time intervals from 1AGL J2022+4032 (left) and 1AGL J2021+3652 (right).The Red line 
indicates the weighted mean of the 42 individual fluxes, from which the^ 2 was calculated: (29 ± 1) x 10" 8 photons cm -2 s" 1 for 
1AGL J2022+4032, and (15 ± 1) x 10" 8 photons cm -2 s" 1 for 1AGL J2021+3652. 



tonic variation in flux with no obvious plateau. The counts maps 
show evidence of the pulsar contribution for all but the lowest 
off-peak exposure fraction, for which the statistical significance 
rfFS =8.4 over two years of data is insufficient to perform a 
variability analysis. For these reasons, a credible analysis of the 
on-peak and off-peak variability based on the AGILE data was 
not feasible. 



4. Discussion of possible counterparts 

In Figure [7] we show the X-ray sources listed in Weisskopf et 



al. 



P006| ) as possible counterparts for 3EG J2020+4017. The 
AGILE contour of the persistent source is consistent with the po- 
sition of the X-ray source [WSC2006] S21, which has been as- 
sociated with the LAT PSR J202 1+4026 w hose pulsations were 
discovered by Fermi. |Trepl et aL] ( |2010| ) searched the XMM- 
Newton archival data and found 2XMM J2021 3 1.0+402645, a 



point source coincident with S21. However, re-analyzing the 
Chandra data, they found that [WSC2006] S25, a strong point 
source within the Fermi 0FGL error box although outside the 
Fermi 1FGL error box, showed evidence of variability during 
the Chandra observation. In addition, S25 is not visible in the 
XMM data, indicating long-term X-ray variability. If S25 is also 
a variable y-ray source, it would be within the AGILE source lo- 
cation accuracy of ~ 1° for the 1-2 week observation durations, 
and could be responsible for the variability apparently observed 
by AGILE below 400 MeV. However, S25 has an infrared coun- 
terpart and could be a normal star. 



The Fermi 11 -month source catalog ([Abdo et al.| 2010a) de- 
tected a nearby source, 1FGL J2020. 0+4049, associated with 
TeV source VER J20 19+407, that is within the one-week error 
circle of AGILE. It was not detected by AGILE, probably be- 
cause of its hard spectrum with index -2.12 + 0.08. Its Fermi 
light curve is consistent with no variability at a level of 33%; 
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Table 2. Variability analysis with and without 10% systematic errors on fluxes. 



Name 


Syst. errors 


x \N df = 41) 


Pvar 


V 


V F 


Pv F 


6F/F 


1AGL J2022+4032, E > 100 MeV 


yes 


66.84 


99.34% 


2.18 


63.26 


0.014 


0.20 




no 


82.45 


99.99% 


3.88 








1AGL J2021+3652, E > 100 MeV 


yes 


40.28 


49.76% 


0.30 


40.03 


0.51 


0.11 




no 


44.47 


67.24% 


0.48 








1AGL J2022+4032, E > 400 MeV 


yes 


48.10 


79.27% 


0.68 


46.07 


0.27 


0.10 




no 


54.46 


92.23% 


1.11 








1AGL J2021+3652, E > 400 MeV 


yes 


34.05 


22.95% 


0.11 


32.77 


0.82 






no 


36.50 


32.94% 


0.17 










Fig. 7. Possible x- ray c ounterparts from 
Weisskopf et al.| ( |2006| ). Green contour 
level: AGILE-GRID 95% C.L. (only statis- 
tical) for 1AGL J2022+4032; cyan circle: 
AGILE-GRID errorbox (statistical + system- 
atic error); red circle: EGRET errorbox for 
3EG J2020+4017; black circle: Fermi errorbox 
for 1FGL J2021.5+4026. 



however, its measured flux for E > 100 MeV decreased from 
(22 ± 6) x 10" 8 photons cm" 2 s" 1 in the first month to below 
4.3 x 10" 8 photons cm -2 s" 1 by the seventh month. An intrin- 
sic flux variability of 34 x 10" 8 photons cm -2 s" 1 over the longer 
AGILE observation period would explain the apparent flux vari- 
ability of 1 AGL J2022+4032. 

Detection of simultaneous flaring in y-ray and, e.g., X-rays 
would provide a definitive identification. However, soft X-ray 
instruments have fields of view too narrow to make triggered ob- 
servations practical. We note that a special class of microquasars 
may not necessarily pr oduce simultaneous hard X-ray emission 
( |Romero& VTTa||2559l hereafter RV09). 

4.1. A variable pulsar 

It is possible that the pulsar LAT PSR J202 1+4026 itself has 
variable y-ray emission for the energy range below 400 MeV. 
This scenario, if confirmed, would open a new field of inves- 
tigation of y-ray pulsars. However, in the absence of other in- 
formation supporting this hypothesis pulsar y-ray variability of 
LAT PSR J202 1+4026 is unlikely. This conclusion is based on 
both previous observational evidence from EGRET and from 
general theoretical consid erations. In particula r, the spin pe- 
riod and time derivative ( |Abdo et all |2009b| ) of LAT PSR 
J2021+4026 CP = 265 ms, P = 54.8 x HT 15 ) assign a unre- 
markable position to this pulsar in both the P - P diagram and 
in the pulsar y-ray luminosity vs. period (or Goldreich- Julian 
current diagram). Instabilities in the y-ray pulsed emission of 
LAT PSR J202 1+4026 might be associated with radio and/or X- 



ray signal changes during the period of the detected y-ray vari- 
ability. Future monitoring of LAT PSR J202 1+4026 can con- 
tribute to test this fascinating hypothesis. 



4.2. Background blazar 

Another possibility is that a blazar behind the Galactic plane 
is contributing variable y-ray emission to the measured flux of 
1AGL J2022+4032. There are a number of known blazars in 
the Cygnus region, but they are well outside the typical AGILE 
one-week error circle. Most optical and radio AGN catalogs 
avoid the Galactic plane because of the high concentration of 
Galacti c sources, heavy extinction, and/or diffuse radio emis- 
sion (|Abdo et aL]|2010b ). Each of these difficulties is especially 



acute at the position of 1AGL J2022+4032, which is located 
within the Gamma -Cygni supernova remnant shell (see Fig. [TJ. 
|Trepl et aL] ( |2010| ) found multiple areas of concentrated radio 
emission within the Fermi error box of 1FGL J2021. 5+4026 as 
well as evidence of variable X-ray emission from nearby source 
[WSC2006] S25. In addition, the y-ray flux could be well be- 
low the detection threshold of AGILE (which, in the presence of 
a known, bright y-ray source is quite high) and still contribute 
to the overall y-ray variability within the one-week AGILE er- 
ror box. Knowing that the intrinsic distribution of blazars should 
be isotropic, we can use t he AGN association s above \b\ > 10° 
in the First AGN Catalog (|Abdo et al.||2010b| hereafter 1LAC) 



to estimate the probability of finding at least one blazar within 
the ^ 1° error box of a one- week observation with AGILE. The 
probability of finding a blazar similar to the 599 1LAC associ- 
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Fig. 5. Simulated cumulative distributions (10000 trials) of the 
variability index V versus fraction of trials, or probability. Solid 
black: mean flux 131 x 10" 8 photons cm -2 s" 1 and square root of 
variance 33.8 x 10" 8 photons cm" 2 s" 1 for E > 100 MeV; dot- 
ted red: mean flux 60 x 10" 8 photons cm -2 s" 1 and square root of 
variance 15.5 xlO -8 photons cm -2 s" 1 for E > 100 MeV; dashed 
green: mean flux 33 x 10" 8 photons cm -2 s" 1 and square root of 
variance 8.5 x 10" 8 photons cm" 2 s" 1 for E > 400 MeV; dot- 
dahsed blue: mean flux 17 x 10" 8 photons cm -2 s" 1 and square 
root of variance 4.4 x 10" 8 photons cm" 2 s" 1 for E > 400 MeV. 
The vertical lines indicate the observed values of V and the as- 
sociated probability for the corresponding dataset (Table [2]). 



ations is « 0.05. A more conservative estimate using only the 
281 1LAC associated FSRQs, yielding a probability of « 0.02 
of finding at least one FSRQ, would better represent the need for 
y-ray variability and the non-detection by Fermi. In either case, 
the probability of chance coincidence is quite low. 

4.3. An X-ray quiet microquasar 

Taking into account both the AGILE-GRID emission above 
100 MeV and the Super- AGILE upper limit in the 15-60 keV 
range (~ 60 mCrab), we consider the possibility that the de- 
tected y-ray variability is caused by transient activity of an X- 
ray quiet microquasar. RV09 analyzed several Galactic sources 
with variable emission in the y-ray energy range and showing 
a ratio L y IL x » 1. They proposed that this kind of emission 
(shortly variable, X-ray quiet, . . . ) can be produced by proton- 
dominated jets in a special class of Galactic microquasars. The 
bulk of the emission at y-ray energies is produced by hadronic 
jets emitted from an accreting source. The model of RV09 pre- 
dicts a y-ray luminosity for this process on the order of « 10 34 
erg/s. Assuming the presence of a y-ray source (an X-ray quiet 
microquasar) within the error box of LAT PSR J202 1+4026, we 
find that that it is required to be at a distance of « 300 pc from 
the Earth, i.e., closer than the pulsar (1-2 kpc). 

The probability of finding this particular type of microquasar 
within the error box of 1AGL J2022+4032 is difficult to quan- 
tify. Nevertheless, because X-ray binaries are concentrated in 
the star-forming regions in the Galactic plane, and high-mass 
X-ray binaries particularly along tangents of spiral arms such as 



the Cygnus region ( [Liu, van Paradijs, & van den Heuvel 2007 
2006), the likelihood that there is an appropriate microquasar 
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Fig. 6. y-ray phase profile of PSR J202 1+4026 as seen by AGILE 
(top) and Fermi (bottom). The AGILE profile is for photons with 
E > 10 MeV while the Fermi profile, taken from |Abdo et al. 
(2009b]), is for E > 300 MeV. 



are isotropically distributed. Similar reasoning applies to such 
po ssible so urce types such as massive s tellar winds (Tav ani et 
[aT| |2009b| ) and novae ( |Abdo et al.||2010d) . 



5. Conclusions 

The apparent y-ray flux variability of 1 AGL J2022+4032 in the 
100-400 MeV range as seen by AGILE leads us to consider two 
possible explanations. One is that the Fermi pulsar, LAT PSR 
J202 1+4026 itself has a variable y-ray flux. This behavior would 
be highly unusual given the properties of y-ray (and radio) pul- 
sars with similar characteristics. Based on the large y-ray un- 
pulsed fracti on of LAT PSR J202 1+4026 as seen by its folded 
light curve ([Abdo et al.| [2009b ), we believe that it is more 



within the error box is much higher than that of blazars, which 



likely that another variable y-ray source within the error box 
of LAT PSR J202 1+4026 contributes to the total emission of 
1 AGL J2022+4032. 1AGL J2022+4032 was not detected in the 
hard X-rays during the AGILE observation periods, with upper 
limits in Super- AGILE varying from 10 to 50-60 mCrab in the 
15-60 keV range. Therefore, if the variable source is a micro- 
quasar rather than an unidentified low-frequency blazar, it would 
be an X-ray-quiet y-ray variety, possibly of the type proposed by 
RV09. 
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We note that during the period in which 1AGL J2022+4032 
was observed by Fermi (starting in August 2008), its flux in 
AGILE-GRID was slightly more stable. In fact, using only the 
AGILE fluxes for E > 100 MeV from the 28 time intervals dur- 
ing the Fermi observations, we find that the variability V was 
1.65 without systematic errors and 0.96 with systematic errors, 
while in the period before Fermi the variability V was 1.94 with- 
out systematic errors and 1.23 with systematic errors. Future 
observations of 1AGL J2022+4032 by both AGILE and Fermi 
will reveal whether this fascinating source continues to show ev- 
idence over the long term. 
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